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ABSTRACT 

Using cosmological simulations with a dynamic range in excess of 10 , we study the transport of 
gas mass and angular momentum through the circumnuclear region of a disk galaxy containing a 
supermassive black hole (SMBH). The simulations follow fueling over relatively quiescent phases of 
the galaxy's evolution (no mergers) and without feedback from active galactic nuclei (AGNs), as part 
of the first stage of using state-of-the-art, high-resolution cosmological simulations to model galaxy 
and black hole co-evolution. We present results from simulations at different redshifts (z = 6, 4, and 
3) and three different black hole masses (3 x 10 7 , 9 x 10 7 , and 3 x 10 8 M Q ; at z = 4), as well as a 
simulation including a prescription that approximates optically thick cooling in the densest regions. 
The interior gas mass throughout the circumnuclear disk shows transient and chaotic behavior as 
a function of time. The Fourier transform of the interior gas mass follows a power law with slope 
— 1 throughout the region, indicating that, in the absence of the effects of galaxy mergers and AGN 
feedback, mass fluctuations are stochastic with no preferred timescale for accretion over the duration 
of each simulation (~ 1-2 Myr). The angular momentum of the gas disk changes direction relative to 
the disk on kiloparscc scales over timescales less than 1 Myr, reflecting the chaotic and transient gas 
dynamics of the circumnuclear region. Infalling clumps of gas, which are driven inward as a result of 
the dynamical state of the circumnuclear disk, may play an important role in determining the spin 
evolution of an SMBH, as has been suggested in stochastic accretion scenarios. 
Subject headings: galaxies: evolution — galaxies: high-redshift — galaxies: nuclei — galaxies: structure 



1. INTRODUCTION 

Supermassive black holes (SMBHs) with masses from 
~ 10 6 to more than ~ 10 9 MfT) are found in the 
centers of most galax i es (e .g., Kormendv fc Richstond 
119951 : iMagorrian et al.l 119981 ) . Observations of correla- 
tions between the masses of SMBHs and several prop- 
erties of their host galaxies indicate that the growth 
of SMBHs is closely tied to th e evolution of their 
hosts (e.g.. IMagorrian et al.l 119981: iFerrarese fc Merritti 
[20001: iGebhardt et al J 120001: iTremaine et al.ll2002D . See 
narios including feedback from active galactic nuclei 
(AGNs) have been explored in order to understand 
both the growth of SMBHs and the nat ure of their 
relati o nship with their host ga l axies ( e.g.. iSilk fc Rees 
l99l| IKauffmann fc Haehne tl 12005 IWvithe fc Loeb 



2001 iDiMatteo et all 120051 iSpringel et all l2005allbl: 



Croton et al.H2006t IDi Matteo et al.H2008l: IHopkins et al l 
20071 Uohansson et al. l l2009bft . AGN s are thought to b e 
powered by accreting SMBHs (e.g.. iLvnden-B cll 1969), 
so that observations of the various phases of AGN ac- 
tivity may provide clues about the connection between 
SMBH g rowth and galaxy evo lution as models explored 
by, e.g., IHopkins et all ((20051) and iSiiacki et all (|2007f ) 
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have shown. It is therefore essential to characterize ac- 
cretion onto SMBHs in order to understand AGNs and 
their relevance for galaxy evolution. 

Substantial SMBH growth requires a mechanism 
for replenishing fuel on the subparsec scales feeding 
the SMBH accretion disk. Large scale tidal torques 
caused by galaxy mergers are an effective mechanism 
for transporting gas from super-galactic and galac- 
tic scales down to scales of several parsecs, where 
other mechanisms may become more viable for fun- 
neling mate rial the rest o f the distance toward th e 
SMBH re.g.. lHernaurs^[l989l:lBa7nes fc Hernauistlll992ft . 
Secular evolution may also play a role as disk in- 
stabilities lead to the formation of bars and spiral 
waves which can transport material dRoberts et al1ll979l : 
iShlosman et~aTl [19891 iRegan fc Teubenl I2004D . There 
are several observations of inflowing gas in the cir- 
cumnuclear regions of low-luminosity AGNs, apparently 
driven by globa l instabilities arising from secular evo- 
lution (see, e.g.. iGarria-Buriilo et al.l 120051: iBoone et al 
2007 1 : iStorchi-Bergmann et all l2007t iLindt-Krieg et al 
2008; Rif fel et alj 12008ft . Mergers are thought to be 
effective for building SMBHs in more massive early- 
type galaxies, which appear to form earlier (see, e.g., 
Hopk ins et al.l I2008 a.b. and references therein), ulti- 
mately building the SMBH-bulge relationships, whereas 
secular evolution may play a larger role in the growth of 
late- type galaxies, particularly at lower redshift (see, e.g., 
IKormendv fc Kennicuttj[2004 and references therein). 

Whichever is the dominant mechanism for fueling 
SMBHs, the specific details of how fueling occurs in 
the circumnuclear regions of galaxies are still not under- 
stood. Fueling may occur continuously, over the course 
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of large-scale dynamical instabilities in the galaxy, or 
it may be an intermittent process, dependent entirely 
on the dynamics on small scales. Cold gas, stochasti- 
cally accreted from the circumnuclear regions of galax- 
ies, ma y be able to sustain the fueli ng of low-luminosity 
AGNs ([Hopkins fc Hernquistl 120061 ). Intermittent ac- 
cretion episodes, consisting of in-falling clouds of gas 
with randomly oriented angular momentum vectors, may 
contribute to the spin-down of SMBHs, thus lowering 
their radi ative efficiency and allowing them to grow 
faster (seelKing fc Pringlel l2006L 120071 : King et al.ll2008f ). 
iWang et al.l (120091 ) have used a Soltan-type argument 
(|Soltanl 11982) to show that the average radiative effi- 
ciency of SMBHs has decreased over time, as might be 
expected if SMBHs grow predominately through episodic 
accretion events that act to lower their spin. Using cos- 
mological simulat ions which fol l owed the buildup of large 
SMBHs at z = 6, iSiiacki et all ()2009f ) showed that black 
hole mergers can contribute to the spin-down of SMBHs 
as well, resulting in low radiative efficiencies and rapid 
growth. A decrease in radiative efficiency must be recon- 
ciled with observations of an increase in th e fraction of 
radio -loud quasars with decreasing redshift (jJiang et al.l 
I2007T ) . which would suggest SMBHs with larger spins 
(and hence larger efficiencies) at low-z. Therefore, un- 
derstanding the behavior of angular momentum in the 
circumnuclear region is an important part of modeling 
accretion, and subsequently AGN feedback. 

Simulations that follow the growth of SMBHs over cos- 
mological times, or over the duration of a galaxy merger, 
often cannot follow the circumnuclear regions of galax- 
ies with high enough resolution to describe the accretion 
flow in detail. Such simulations must make approxima- 
tions for the accretion rates using the properties of the 
galaxies on scales that are resolved. A common technique 
is to assume t hat the unresolved disk is fed by Bondi- 
type accretion (IBovle fc Lvttletonlll939t iBondi fc Hovld 
1944; Bondi 1952), as in for exam ple the smoothed parti- 
cle hydrodynamics simulations of Spring eTet al.l ([2005a) 
who estimate the accretion rate based on the properties 
of the gas on a scale of ~ 100 pc. The assumption of 
Bondi accretion appears to be reasonable for following 
evolution over cosmological times, where the average ac- 
cretion rate onto the black hole cannot be more than 
the average accretion rate through any radius, since the 
fuel supply is limited by large scales. Previous cosmo- 
logical simulations have been successful in reproducing 
observed population demographics and trends, such as 
the black hole mass f unction and galaxy colors and mor- 
phologies (as in, e.g., Pi Matteo et al.l (20081 : iCroft et ail 
120091: (McCarthy et al.ll2009l) . However, for detailed stud- 
ies of the growth and evolution of individual SMBHs or of 
the coupling between AGN feedback and the black hole 
accretion rate, more accurate descriptions of the accre- 
tion rate and its dependence on the small-scale features 
of the host galaxy are desirable. 

Small-scale simulations have addressed gas dynam- 
i cs in subgalactic-scale disks with high resolution 
Fukuda et aLll2000l IWadalf200ll: iWada fc Normanll2001l: 



Escalall2007l : IWada fc Normanl 120071 : iKawakatu fe Wadal 
20081) . finding the development of a turbulent, multi- 



phase interstellar medium. The approximation of Bondi 
accret ion in such an environ ment is not necessarily in- 
valid. iKrumholz et al.l (|2005l) have shown that modified 



forms of the Bondi prescription can describe accretion in 
turbulent environments. A simulation must be equipped 
to model the properties of the turbulence in order to 
employ the modified Bondi prescription. If the gas con- 
tains a significant amount of angular momentum, the un- 
modified Bondi prescription gives inaccurate estimates 
of the accretion rate as well (jProga fc Begelmanl 120031 : 
IKrumholz et al.l I2006D . However, even modified Bondi 
prescriptions become inapplicable in the case of the self- 
gravitating, rotationally supported disks that are likely 
to form as large amounts of gas are driven inward in 
high-redshift galaxies. As the use of adaptive techniques 
increasingly improves the resolution of cosmological sim- 
ulations, accretion onto black holes can no longer be de- 
scribed by approximate prescriptions such as Bondi ac- 
cretion. 

In the present paper, we use cosmological adaptive 
mesh refinement simulations with a large dynamic range 
to study the transport of gas and angular momentum 
through the circumnuclear disk of an SMBH host galaxy 
over time. The goal is to provide a description of ac- 
cretion that can be compared to prescriptions typically 
applied in larger scale simulations. It will be shown that 
in the limiting case of relative quiescence (e.g., in be- 
tween merger events and without AGN feedback) accre- 
tion is a stochastic rather than a continuous process. The 
method behind the simulatio ns used here is described in 
detail in iLevine et al.1 (|2008f ). hereafter Paper I, and in 
iLevind (|2008| ). In Section we briefly summarize the 
details of the simulations. The results of the simulations 
are given in Sections [5][5l including an analysis of the 
mass accretion rate and the angular momentum of the 
gas in the circumnuclear region of the galaxy. Finally, 
we summarize and discuss the results in Section [6l 

2. SIMULATION 

The simulations presented here were run using the 
Adap tive Refinemen t Tree (ART) code dKravtsov et all 
ll997tlKravtsovlll999l:lKravtsov et al.ll2002D . following the 
"zoom-in" method described in detail in Paper I. The 
code follows gas hydrodynamics on an adaptive mesh 
and includes dark matter and stellar particles. The gas 
cooling and heating rates are tabulated as functions of 
density, temperature, metallicity, and redshift over the 
temperature range 1 2 < T < 10 9 K using CLOUDY 
(jFerland et al.1 119981 ) , which accounts for the metallic- 
ity of the gas, the formation of molecular hydrogen and 
cosmic dust, and UV heating due to cosmological ioniz- 
ing background. Because the circumnuclear disk of the 
simulated galaxy reaches high densities, the gas can be- 
come optically thick to its own cooling radiation. In or- 
der to account for this effect and approximate optically 
thick cooling, we have also run a simulation that uses a 
Sobolev-like approximation for the column density of the 
gas to account for the large opacity of the highest-density 
regions in the simulated galax y in calculating the cooling 
rates (see IGnedin et al1l2009l) . 

Stellar particles are formed in cells with appropriate 
densities and temperatures at an efficiency that matches 
observed star f ormation rates on kilop arsec scales 
(iKennicuttl 119981) and 10 pc scales (e.g., lYoung et al.l 
19961 iWong &"Bl itz 2002). Presently, only the cosmo- 
logical portion of the simulation (before zooming-in) in- 
cludes radiative transfer and feedback and enrichment 
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TABLE 1 

Summary of Each Simulation Run 



Run 


2 


A/200 




M 4 


M. 


Ax 


Duration 






(lO 11 /!" 1 M ) 


(lO^fe^Ms) 


(lO 9 /!" 1 M ) 


(1O 7 M ) 


(proper pc) 


(10 6 yr) 


Z3L20 


3 


2.4 


3.3 


7.3 


3 


0.032 


0.84 


Z3.5L20 


3.5 


2.0 


2.8 


5.3 


3 


0.028 


0.71 


Z4L15 


4 


1.7 


2.5 


4.0 


3 


0.819 


1.0 


Z4L18 


4 


1.7 


2.5 


4.0 


3 


0.102 


1.0 


Z4L20 b 


4 


1.7 


2.5 


4.0 


3 


0.026 


1.9 


Z4L20.B3 


4 


1.7 


2.5 


4.0 


9 


0.026 


1.0 


Z4L20.B10 


4 


1.7 


2.5 


4.0 


30 


0.026 


0.80 


Z4L20.OT c 


4 


1.7 


2.5 


4.0 


3 


0.026 


0.85 


Z6L20 


6 


0.66 


1.0 


0.84 


3 


0.018 


1.9 



a Mass of gas with T < 20,000 K. 
b Fiducial run. 

c Includes a Sobolev-like approximation for the column density to model optically thick cooling. 



from stars (although stellar particles are still formed in 
the zoom-in simulations). The physics on small scales 
in galaxies is complex, so our approach is to individu- 
ally address each problem, carefully building a realistic, 
state-of-the-art simulation one piece at a time. 

Here, we provide a brief overview of the zoom-in 
method, with a complete description given in Paper I. 
A cosmological simulation is evolved from a random re- 
alization of a Gaussian density field at z — 50 in a 
periodic box of 6 /i -1 Mpc with an appropriate power 
spectrum, and is followed assuming a flat ACDM model: 
O = 1 - = 0.3, n h = 0.043, h = #o/100 = 0.7, 
n s = 1, and er 8 = 0.9. Beginning with the cosmo- 
logical simulation, which has a maximum resolution of 
w 180/i _1 pc (37/i _1 pc proper at z = 4), the resolution 
is slowly increased one refinement level at a time, reach- 
ing a quasi-stationary state on each level before increas- 
ing to the next level. The final, maximum resolution in 
our fiducial run is ~ 0.089/i _1 pc (0.02/i~ 1 pc proper at 
z = 4), corresponding to 20 levels of refinement. 

After reaching the maximum resolution, a fraction of 
the gas in the center of the galaxy is replaced with a black 
hole particle of equal mass and momentum. The circum- 
nuclear region is extremely gas rich, with the gas com- 
prising 91%, 83%, and 79% of the baryon mass inside a 
radius of 1 kpc at z = 6, 4, and 3, respectively. In this re- 
gion, the gas mass dominates over all other components, 
including the black hole particle. The black hole mass 
has not been chosen to match any of the local scaling re- 
lations (e.g.. iMagorrian et alj|1998t iFerrarese fc Merritti 
12005 iGebhardt et all 120001 : iTremaine et al.l 120021 ). In- 
stead, we have chosen a mass that is large enough to dis- 
tinguish it from stellar particles in the simulation, and 
that tends to stay put in the center of the galaxy, but 
does not depend on any a priori assumptions about how 
the scaling relations extend to high redshift. The pos- 
sibility remains that the simulated galaxy should host a 
more massive black hole, but since we do not follow the 
black hole's growth from its seed formation, the mass of 
the black hole particle is somewhat arbitrarily chosen. 
Presently, the mass of the black hole particle does not 
change over the course of the zoom-in simulations. This 
simplification allows us to avoid making any assumptions 
about how the gas mass on the smallest scales of the sim- 
ulations will be accreted onto the black hole particle. 

After the introduction of the black hole particle, the 



simulation continues to evolve with the maximum reso- 
lution for several hundred thousand years. The duration 
of each simulation varies with redshift. At lower red- 
shifts, the increased mass of the galaxy results in more 
high-density cells, requiring high resolution. Following 
a larger portion of the simulation with finer resolution 
is more computationally expensive, therefore the lower 
redshift simulations take longer to evolve. As a result, 
e.g., the z = 3 zoom-in is evolved only half as long in 
physical time as the fiducial z = 4 simulation. 

Most of the following analysis is applied to the short- 
term evolution of the galaxy (~ 10 6 yr) during a zoom-in 
episode at z — 4 containing a 3 x 10 7 M Q SMBH parti- 
cle (considered the fiducial ruifl). Additional simulations 
were performed, including zoom-in episodes at other red- 
shifts with the same black hole mass (z — 3, 3.5, and 
6), runs with different black hole masses (9 x 10 7 and 
3 x 10 s M©) at z = 4, two lower resolution versions of 
the fiducial run (with maximum levels of refinement 15 
and 18), and a version of the fiducial run that implements 
an approximation for optically thick cooling. The z — 4 
simulation was chosen for the fiducial run because this 
redshift corresponds to a relatively quiescent stage in the 
growth of the host galaxy. At this stage, the main galaxy 
of the simulation has not undergone a large mass ratio 
merger in ss 400 Myr, although the galaxy is still actively 
growing. At z = 6 (also studied here), the galaxy is less 
relaxed, as it is entering into a 4:1 merger at that time. 
It is expected that the mass of the black hole should con- 
tinue to grow with decreasing redshift. Since the differ- 
ent redshift simulations each contain the same black hole 
particle mass, they do not provide the growth history 
of a single galaxy, but rather they represent the distinct 
growth histories of similar disk galaxies. Table [T] summa- 
rizes the details of the simulations performed, including 
redshift, the total mass within a radius enclosing a mean 
density of 200 times the present-day critical density of 
the universe, centered on the main galaxy (M200), the 
mass of gas with temperature less than 20,000 K (M g ), 
the stellar mass (M*), the black hole mass (M,), the size 
of the most refined cell in each run (Aa; m i n ), and the 
duration of each run in physical time. 

7 The fiducial run, Z4L20, is effectively the same simulation 
described in Paper I. The only difference is in the initial zoom-in of 
the run presented here, which occurred more slowly (see discussion 
in Section 2.2, Paper I). 
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Fig. 1. — Time-averaged interior gas mass as a function of radius 
for each simulation. 

Paper I gives the details of the spatial structure of 
the circumnuclear region for the fiducial run. Each run 
described here develops a similar structure: a cold, rota- 
tionally supported, self-gravitating gas disk develops in 
the circumnuclear region of the galaxy (inside ~ 100 pc), 
with a steep power-law density distribution. The disk is 
globally unstable, leading to the development of waves 
and instabilities, which drive turbulence on a range of 
scales, injecting energy through shocks. The turbulence 
is highly supersonic, effectively raising Q and acting to 
maintain local stability in the disk. In the run Z4L20.OT, 
the approximation for optically thick cooling accounts 
for the densest gas being optically thick to its own cool- 
ing radiation. Consequently, the mass-weighted mean 
temperature in that run is higher than in the fiducial 
run, which contributes to a slightly higher Toomre Q- 
parameter. In cither case, the higher effective Q prevents 
the circumnuclear disk from fragmenting entirely into 
star-forming clumps on a free-fall time. However, the ac- 
tual star formation rate in the simulation does not reflect 
the complicated behavior of the circumnuclear region, 
because the rate is simply density- and temperature- 
dependent and operates on timescales much longer than 
the million-year time scales of the zoom-in simulation. 
During the zoom-in simulation, the disk shows transient 
features caused by spiral waves and global instabilities, 
which slowly allow angular momentum transport, driv- 
ing gas inward toward the center of the galaxy. Sim- 
ilar behavior, where global instabilities generate stabi- 
lizing turbul ence, has been observed i n other numerical 
simulations (iRegan fc Haehneltl [200l iWise et "ail 120081: 
lEnglmaier fc Shlosman 20041) and h a s been studied an- 
alvticallv by iBegelman fc Shlosmanl f|2009T ) . In the fol- 
lowing sections, we explore the behavior over time of the 
mass accretion rate and the angular momentum of the 
gas throughout the circumnuclear disk. 

3. DISK MASS AND THE RATE OF TRANSPORT 

Figure [1] shows the mean gas mass, M (averaged over 
the duration of the simulation), interior to radius, r, for 
each simulation (except z — 3.5). The radius r is mea- 
sured with respect to the position of the black hole par- 
ticle, which approximates the center of the galaxy over 
the duration of each zoom-in simulation (except for a few 
brief intervals, which are discussed later in this section). 



A comparison of the average gas mass in each zoom-in 
episode approximately shows the growth of the circum- 
nuclear disk over cosmological times. The substantial 
growth in mass between z = 6 and z — 4 is caused by a 
major merger (with a mass ratio of 4:1) and a few minor 
mergers occurring in between. 

In contrast with the averaged mass-profiles, the in- 
stantaneous mass-profiles at individual times during the 
zoom-in part of the simulation (not shown) reveal the 
transient nature of the circumnuclear disk. As the glob- 
ally unstable disk forms transient structures on a range of 
scales, the interior gas mass can vary significantly. Fig- 
ure[5]shows the amplitude of mass fluctuations, M(t)/M 
at four different radii (ranging from 0.1 to 100 pc, span- 
ning scales which have undergone at least one rotation 
period over the course of each of the zoom- in simulations) 
over time for each of the three redshift simulations. The 
figure gives an indication of the turbulent nature of the 
gas in the circumnuclear disk of the galaxy. In their 
high-res olution re-simulations of m ergers between disk 
galaxies. iHopkins &; Quataer t (2009) find a qualitatively 
similar picture, in which global instabilities operate on a 
range of scales, causing a great deal of fluctuation in the 
interior gas mass as a function of time throughout the 
circumnuclear disk. 

Around 550,000 yr in the z — 4 simulation, a ~ few 
xl0 7 M Q clump of gas (which forms in the central 10 pc, 
perhaps as a result of the global instabilities in the disk) 
moves within a few parsecs of the black hole particle, 
temporarily creating a gravitational disturbance that dis- 
places the black hole particle for w 50,000 yr. The black 
hole soon resettles into the bottom of the potential well 
of the disk, but because the interior mass is measured 
with respect to the black hole particle, the interior gas 
mass measurements are temporarily affected by the black 
hole's motion (see the dotted lines in Figure [2]). There- 
fore, the interior gas mass shown in Figure [5] has been 
interpolated over the duration of the black hole displace- 
ment in all subsequent computations (unless otherwise 
noted). The properties of the clump of gas may be sig- 
nificant for the angular momentum of the black hole ac- 
cretion disk, which will be discussed further in Section 
El Runs Z4L20.BH3, Z4L20.BH10, and Z4L20.OT un- 
dergo similar events of variable duration. The details 
of each displacement, such as the size and the duration 
of the subsequent relaxation, depend on the somewhat 
arbitrarily chosen black hole particle mass and the gas 
dynamics, which vary between different runs. Therefore, 
even though the effect of in-falling clumps on the gas disk 
may be significant, the displacement of the black hole 
particle in these simulations does not necessarily imply 
that such displacements occur similarly in nature. 

As Figure [2] shows, a characteristic accretion rate is 
not straightforwardly determined by any one individual 
snapshot of the circumnuclear region, given the amount 
of fluctuation in the gas mass over the timescale of the 
simulations. A simple estimate of the net accretion rate 
could be obtained by measuring the difference in mass 
at each scale between the beginning and the end of each 
run. However, the dynamical and turbulent behavior 
of the circumnuclear disk gives rise to a highly variable 
accretion rate, often indicating as much outflow as inflow 
of gas mass, and in some instances, resulting in a net 
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0.0 0.5 1.0 1.5 0.0 0.5 1.0 1.5 
t [ 1 6 yr] t [1 6 yr] 

Fig. 2. — Amplitude of gas mass fluctuations at redshifts 3, 4, and 6 for four different radii. The dotted black line shows the fluctuations 
measured with respect to the black hole particle position over the duration of the black hole displacement. The thin solid line shows 
M = M for comparison. (Figure [l] shows M as a function of radius, r, for each simulation). 

outflow. In this Fourier transform of the interior 

gas mass provides a natural characterization of the mass 
fluctuations and their dependence (if any) on timescale. 

We have chosen the sampling frequency for the mass 
(and all other physical quantities), somewhat arbitrar- 
ily, to correspond to one time step on level 6. The fre- 
quency is such that the entire zoom-in region in each of 
the simulations steps at least 8 times (corresponding to 
the least-resolved level of refinement in the circumnuclear 
disk; level 9), and at most 16,384 times (corresponding 
to the most-resolved level; level 200- During the sim- 
ulation, the physical time corresponding to each step is 
redetermined after every level-0 step, and depends on the 
dynamics of the disk (so that the time step obeys a global 
CFL condition). Therefore, the length of time between 
samplings varies throughout the duration of the simu- 
lation. The Lomb-no rmalized periodogram (see iLombl 
ll976HPress et al.lll992l Numerical Recipes, Section 13.8), 
or the spectral power normalized by the variance of a 
data set, is useful for studying the behavior of unevenly 
sampled data, such as the uneven time steps here. The 
normalized spectral power as a function of the angular 
frequency, uj, is given by 



-i 2 



XX-Wj - M) cos u>(tj - t) 
i 



2°m cos2 — r ) 



+ 



- M) smui(tj - t) 

j 



2cr M J2 s ' m - T ) 

j 



(1) 



summed over N data points, where cr^j is the variance in 
the interior gas mass. The quantity r, defined by 



tan(2wr) 



N 

sin 2utj 
j=i 

N 

cos2ujtj 

i=i 



(2) 



8 Time steps forward in the simulation so that for each step on 
level 6, the simulation has taken 2' -6 steps on level I. 



is introduced so the spectral power is independent of 
shifting the data in time. Although the Lomb-normalized 
periodogram is not directly computed from the Fourier 
transform of the data, we use it to estimate the Fourier 
transform of the interior gas mass as follows: 
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Fig. 3. — Mass flux, estimated from the spectral power of the interior gas mass fluctuations, for three different maximum resolution 
simulations (f5, f8, and 20 levels of refinement) at z = 4. Here, the flux is only calculated for the first million years of run Z4L20, for a 
consistent comparison with runs Z4L15 and Z4L18. 
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The Fourier transform of the mass follows a power law 
with slope « — 1 in each of the simulations, showing no 
substantial departure down to frequencies comparable to 
the rotational period of the disk. The time derivative of 
the mass Fourier transform gives a characteristic "accre- 
tion rate" through radius r as a function of frequency, 



M(u) = iwM(w). 



(4) 



Strictly speaking, the accretion rate given by Equation 
H]does not describe the net inward transport of gas, but 
rather it gives typical values of the flux of gas mass into 
and out of a given region on a timescale 2tt/lo. The ab- 
solute value of the flux given by Equation |4] is shown 
as a function of frequency in Figures [3][5] The data are 
smoothed over neighboring points (the number of neigh- 
bors used in the smoothing is proportional to angular fre- 
quency, ui) to reduce the noise so that the different curves 
are readily distinguishable. Figure |3] shows a comparison 
of the z = 4 runs with different maximum resolution 
(Z4L15, Z4L18, and Z4L20). The mass flux in each run 
is only shown for scales that are resolved by four cells or 
more. We find good agreement in the behavior shown 
for different resolution simulations, indicating that the 
results are converged with respect to resolution. 



The mass flux is larger for decreasing redshift, because 
the galaxy itself is more massive, but the flux follows the 
same stochastic behavior at different epochs (Figur^lJ), 
as well as for different black hole masses and for the op- 
tically thick cooling run (Figure [5]). We note that 0.1 pc 
(proper) is the resolution limit at z = 4 (corresponding 
to four level-20 cells), and therefore this scale is barely 
resolved at z = 4 and under-resolved at z = 3. There- 
fore, the top-left panel of Figure [4] does not show the 
z = 3 simulation. The vertical lines in each panel corre- 
spond to the average angular speed of the gas at a given 
scale in each simulation. The fluctuations do not appear 
to be correlated with the orbital period at any radius in 
the circumnuclear disk. The flat slope in the flux, M(ui), 
shows that the rate of transport is independent of fre- 
quency, i.e., that accretion is a stochastic process with 
no preferred timescale. The result is generic over each 
run studied here, including the optically thick cooling 
run, which is dynamically similar to the other runs de- 
spite the presence of higher temperature gas throughout 
the circumnuclear disk. Simulations of self-gravitating 
disks show behavior similar to our own results, with 
instabilities driving turbulence (Wada fc Normanl 119991 : 
iWadal I2001t iWada et all 120021 ), which ultimately leads 
to intermitte nt fueling episodes varying on ti mescales of 
10 4 - 10 5 yr (IFukuda et allllM IWaHall200l . 



4. COMPARISON TO OTHER ACCRETION 
PRESCRIPTIONS 

As mentioned in the Introduction, simulations that do 
not resolve the accretion disk of an SMBH, or the re- 
gion where the gravity of the SMBH dominates, typically 
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Fig. 4. — Mass flux, estimated from the spectral power of the interior gas mass fluctuations, for three different rcdshifts. The vertical 
lines correspond to the mean angular speed of the gas disk at radius, r, in each simulation. 
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Fig. 5. — Mass flux, estimated from the spectral power of the interior gas mass fluctuations, for three different black hole masses and 
for the optically thick cooling run at z = 4. The vertical lines correspond to the mean angular speed of the gas disk at radius, r, in each 
simulation. 



include some prescription for the accretion rate that de- 
pends on the properties of the simulated galaxy on scales 
that are resolved. It is useful to check how well a prescrip- 
tion, such as Bondi accretion, describes mass transport 
when extended down to parsec scale resolution. 

The original Bondi prescription assumes pressure- 
supported gas in hydrodynamic equilibrium. Turbu- 
lence, as well as angular momentum (such as in a 
disk), has been shown to lower the actual accretion 
rate below the Bondi rate (lAbramowicz fc Zurekl 119811 : 
iProga fc Begelmanl 120031: iKrumholz et all I2006T ). The 
circumnuclear disk in the simulations here is both turbu- 
lent and rotationally supported (see Figure 6 of Paper I), 
bearing little resemblance to the medium modeled in the 
original Bondi prescription. Modifications to the Bondi 
prescription that account for both turbulence and vortic- 
ity in gas can be applied to estimate t he reduced accre- 
tion rate (jKrumholz et all [2001 [20061 ). However, even 
such a modified prescription fails to describe accretion of 
self-gravitating gas, such as in the circumnuclear disk of 
the present simulations. Disks that are self-gravitating 
form bars and spiral waves that further change the flow 
of gas. 

Figure [5] shows a comparison of the Bondi Hoyle rate 
for turbulent gas (computed using the spherically av- 
eraged radial profiles of the gas density, turbulent ve- 
locity, and sound speed) and the vorticity-reduced rate, 
applied to each redshift simulation. In calculating the 
Bondi Hoyle rate for turbulent gas, A/bh, the thermal 
sound speed of the gas has been replaced by an ef- 
fective sound speed defined as the quadrature sum of 
the thermal sound speed and turbulent velocity. The 
vorticity-reduced rate. M tU rb, given by Equation 3 of 
IKrumholz et al.l (|20 06^. is computed cell by cell, pro- 
ducing a log-normal probability density function (PDF) 
of accretion rates within individual spherical shells cen- 
tered on the black hole particle. The geometric mean of 
the distribution, (M tur b), gives characteristic accretion 
rates under this prescription and is shown for several 
different radii in Figure El The vorticity-reduced rates 
are 2 — 3 orders of magnitude smaller than the standard 
Bondi Hoyle rates. These results are qu alitatively con- 
sistent with those of iDeBuhr et a l. ( 2009) , who find that 
accretion in their merger simulations is regulated by an- 
gular momentum transport processes in the host galaxy, 
keeping the accretion rate below the Bondi rate, except 
during peak activity (when the accretion rate approaches 

M E dd, which is M B h)- 

As revealed in Figure [6l the Bondi prescription gives 
large estimates for the accretion rates (3> lOAiEdd) inside 
the central few hundred parsecs. The steep, power-law, 
density profile of the gas in the circumnuclear disk con- 
tributes to the large accretion rates shown in Figure [6] 
(which are proportional to density). The inclusion of an 
approximation for optically thick cooling in the simula- 
tions does not significantly change the density profile or 
the velocity dispersion of the gas. Therefore, the Bondi 
prescription produces similar results for runs Z4L20 and 
Z4L20.OT (dashed and solid curves in Figure [6]). Ad- 
ditional effects not included in our current simulations 
may contribute to the depletion of gas in the circumnu- 
clear disk, such as AGN feedback and stellar feedback 
(not included in the zoom-in portion of the simulation), 
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Fig. 6. — Bondi accretion rate as a function of radius for dif- 
ferent redshifts. Two methods of estimation are shown: (1) Mbh 
computed from spherically averaged quantities (solid black curve: 
fiducial run, dashed curve: optically thick cooling run) and (2) the 
geometric mean and standard deviation of the M tur b PDF (points). 
The Eddington accretion rate, MEdd f° r a 30 million Mq black 
hole with a radiative efficiency of 0.1 is also shown for comparison 
(dotted black line). 



thus lowering the estimated Bondi rate. However, if the 
disk remains self-gravitating and therefore susceptible to 
instabilities, the Bondi prescription will inaccurately de- 
scribe the transport of gas through this region. 

The average accretion rate through the circumnuclear 
disk over cosmological times can be estimated by com- 
paring the mean interior gas mass (the averages shown 
in Figure [lj for the different redshift simulations. Fig- 
ure [7] shows the mean gas mass, interior to radius r as 
a function of the age of the universe, t agc . Once again 
we emphasize that the different redshift simulations do 
not necessarily describe different stages of growth of the 
same galaxy because they each contain the same mass 
SMBH (rather than a black hole that grows with red- 
shift). However, the black hole particle does not cur- 
rently play a large role in the evolution of the simulated 
galaxy (at least not below z = 4, where the mass of the 
black hole is dominated by the gas mass all the way down 
to the resolution limit). 

The dashed line in Figure [7] shows the mass of a black 
hole, initially 3 x 10 7 Mq at z = 6, if it grows continuously 
at the Eddington limit according to 



M(t) = M, e 



t/t s 



(5) 

where t$ is the Salpeter time (jSalpeterl Il964h of 4.5 x 
10 7 yr, for a radiative efficiency r\ = 0.1. If the dynam- 
ics of the circumnuclear disk extend all the way down to 
scales beneath the resolution (which cannot be assumed) 
then the growth rate approximated by the simulation 
points, which appears to steepen at early times with de- 
creasing scale, may continue to steepen down to smaller 
scales, approaching the Eddington limit. Therefore, a 
black hole in our simulations will be able to grow accord- 
ing to Equation [5] at high long as efficient fueling 
from the circumnuclear disk can be sustained. This is 
perhap s consistent with the is olated disk galaxy simula- 
tions of iSpringel et alj (|2005af ). where a sufficiently mas- 
sive black hole will grow according to Equation [5l in the 
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Fig. 7. — Mean gas mass interior to radius r as a function of the 
age of the universe, t age . The error bars correspond to the disper- 
sion in the interior gas mass over each run. A z = 3.5 simulation, 
not mentioned above, is also included in the figure. 

absence of AGN feedback. Once gas reaches sub-parsec 
scales, the actual accretion rate onto the black hole is 
governed by the physics of the accretion disk, not mod- 
eled by our simulations. Nonetheless, Figure [7] shows 
that the increase in the mass of the circumnuclear gas 
disk over cosmological times yields enough fuel to build 
up the mass of the SMBH if the fuel can be efficiently 
accreted. 

5. ANGULAR MOMENTUM 

The large fluctuations in the interior gas mass corre- 
spond to waves moving through the circumnuclear disk, 
which can potentially drive individual clumps of molec- 
ular gas into the vicinity of the SMBH. The fluctuations 
develop stochastically, as shown in Section [3J which may 
ultimately lead to stochastic SMBH accretion events. If 
the direction of the angular momentum vector of ac- 
creted material with respect to the black hole's spin 
axis also varies stochastically, some fraction of accretion 
events will be anti-aligned, lowering the black hole's spin 
and consequently its radiative efficiency (|King fc Pringld 
2006, 2007$ INavakshin fcKmel l20"07t Km7 et all I2008D . 
However, large amounts of accretion, as in merger driven 
fueling, will tend to align an SMBH's spin with that of 
the circumnuclear disk , increasing its radiative efficiency 
(jVolonteri et al.lfeOOTi ). Spin-down resulting from accre- 
tion could be more effective in non-merging disk galaxies 
hosting smaller SMBHs, as secular processes stochasti- 
cally drive small amounts of gas toward the SMBH. The 
resulting low spin may be consistent with observations 
of disk galaxies as hosts for radio-quiet sources, if radio- 
loudness is associated w ith the spin of the black hole 
(|Berti fc Volonteri| [2008). Since black hole mergers may 
also c ontribute to the spin-down of SMBHs (|Siiacki et al.l 
2009), it is important to qualify the relative importance 
of gas accretion for determining black hole spin. 

The simulations here do not model the gas down to 
scales inside the accretion disk, nor allow the black hole 
particle to accrete material from larger scales (as dis- 
cussed in Section [2]). For these reasons, the simulations 
do not follow the spin of the black hole. The simulations 
do model gas on parsec scales, enabling us to resolve the 



direction of the angular momentum vector of gas and to 
provide boundary conditions for simulations of SMBH 
accretion disks. 

The angular momentum vector is measured with re- 
spect to the black hole particle's position and velocity. 
Some uncertainty in the direction of the angular momen- 
tum vector could arise because the position of the black 
hole does not necessarily coincide exactly with the po- 
sition of the center of the galaxy and because the black 
hole may have a significant rotational velocity about the 
center. During most of the simulation, the black hole 
particle lies in the densest cell, which coincides with the 
minimum of the galaxy's gravitational potential. There- 
fore, the black hole particle typically provides a reason- 
able estimate of the location of center of the galaxy. We 
have also measured the angular momentum relative to 
a center whose velocity is defined by an average over 
the central cell and six of its neighbors, but find that 
the effect on the angular momentum vector is negligible 
(except for the brief periods of black hole particle dis- 
placement). 

Figure [5] shows a map projection of the direction of 
the normalized angular momentum vector as it evolves 
in time at 100, 10, and 1 pc from the black hole particle 
for the z = 4 and z = 3 simulations (left and right, re- 
spectively). The mean direction of the rotation axis of 
the galaxy on kiloparsec scales, which remains compara- 
tively constant over the ~ 1 Myr duration of the zoom-in 
simulation, is oriented toward the center of the map in 
each case. The disk is slightly warped on scales smaller 
than a kiloparsec, so that the axis of the circumnuclear 
disk is oriented at an angle to the large-scale disk. Figure 
[9] shows the inclination angle of the angular momentum 
vectors shown in Figure [8] with respect to those mea- 
sured at a kiloparsec, over time. Together, Figures [8] 
and [9] demonstrate the effect of the circumnuclear disk's 
chaotic behavior on the angular momentum of the ac- 
creting gas. 

At z = 4 and a scale of 100 pc, the direction of the 
angular momentum vector changes slowly but shows lit- 
tle scatter over the course of the simulation. There is a 
small amount of scatter during the displacement of the 
black hole particle, which is visible on each scale shown. 
Rather than interpolating the magnitude and direction 
of the angular momentum vector of the gas over the du- 
ration of the black hole displacement, as was done with 
the interior gas mass profile in Figure [21 we continue 
to measure quantities relative to the black hole particle. 
Once the black hole particle resettles into the bottom of 
the potential well, the excess scatter in L/L vanishes. 
The slow change in direction at 100 pc may correspond 
to the increased warping of the disk as the simulation 
progresses. In a similar an alysis of simula t ions o f the 
growth of a disk galaxy by iSaitoh fc Wadal (|2004f >. the 
spin axis of the central gas disk also changed direction 
on 100 pc scales (relative to kiloparsec scales). 

The angular momentum vector at 10 and 1 pc starts 
out with the same approximate orientation as the larger 
scales but shows slightly more scatter. There are two 
sudden changes in the direction of the angular momen- 
tum vector at the 1 pc scale. In the first incident, at 
~ 0.55 Myr, the angular momentum changes abruptly 
by more than 100°, coinciding with the time of the dis- 
placement of the black hole particle described in Section 




Fig. 8. — Map showing the direction of the normalized angular momentum vector, L/L as it changes over time at 100, 10, and lpc in the 
2 = 4 and 2 = 3 simulations. The colors correspond to the time elapsed in the simulation since the introduction of the black hole particle. 
Each map is aligned so that the mean rotation axis of the galactic disk at a scale of lkpc lies in the center. This figure is best viewed in 
color. 




Fig. 9. — Inclination angle of the angular momentum of the inner disk with respect to that of the kiloparsec disk over physical time, for 
2 = 4 and 2 = 3 (left and right panels, respectively). 



[3] (as a clump of gas forms near and moves into the cen- 
ter). The flip in direction is also visible at 10 pc. In the 
second incident, at ~ 1.2 Myr, the angular momentum 
shifts by ~ 100° while the black hole remains stationary 
at the center of the galaxy. The second shift is the result 
of gravitational interaction with a ~ few x 10 8 M Q clump 
of gas, which develops at > 100 pc and continues to move 
toward the center of the disk at late times. The mass of 



the clump is comparable to the mass of the disk interior 
to ~ 10 pc. The second shift is also visible at the 10 pc 
scale but occurs less suddenly. The two panels in Fig- 
ure [10] show three-dimensional volume-renderings of the 
gas density at 30 pc before and after the flip (top and 
bottom, respectively). The inset shows a zoomed-out 
view of the disk, where the large clump of gas is visible 
> 200 pc from the center. If the clump remains intact 
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Fig. 10. — 3D volume-rendering of the gas density before (top) 
and after (bottom) the second angular momentum shift in the fidu- 
cial run referred to in the text (and in Figures 151 and l9)l . The inset 
in the top panel shows the large clump developing at a few hundred 
parsecs that is influencing the dynamics of the inner circumnuclear 
disk. The circle in the center of each panel corresponds to the 
position of the black hole particle. This figure is best viewed in 
color. 

as it reaches the center of the circumnuclear region, a 
corresponding temporary displacement of the black hole 
particle is expected, similar to the earlier displacement. 
The black hole particle is likely to be displaced when 
the in-falling gas mass is comparable to or greater than 
the black hole mass and the mass of the gas disk inte- 
rior to the in-falling clump. However, a displacement has 
not yet occurred by the end of the simulation. Similar 
angular momentum flips occur over the course of runs 
Z4L20.BH3, Z4L20.BH10, and Z4L20.OT (not shown), 
often corresponding to coincident displacements of the 
black hole particle. 

The z — 3 simulation shows more scatter at 100 pc than 
the z = 4 simulation. There is a rapid shift in the direc- 
tion of the angular momentum vector of ~ 30 — 45° be- 
ginning at ~ 0.35 Myr in the z — 3 simulation. The shift 
is less dramatic than both of those observed in the z = 4 
simulation. There is no corresponding black hole particle 
displacement, likely because the mass of in-falling clump 
is small compared to the mass of the disk (and the black 
hole). The orientation of the disk appears to be more 
coherent between the different scales in z — 3 simulation 
than in the z = 4 simulation. This may be a result of the 



circumnuclear disk being more massive at lower redshift. 
In broad terms however, the behavior of the angular mo- 
mentum vector in each of the simulations consistently 
shows that the gas delivered to the SMBH may have 
varying angular momentum, which can shift suddenly 
as massive clumps of gas develop and move through the 
disk. The properties of the gas that is ultimately accreted 
may then determine the spin and subsequently the mode 
of feedback produced by the SMBH. 

6. DISCUSSION AND CONCLUSIONS 

We have studied the transport of gas in the circumnu- 
clear region of a disk galaxy within a cosmological simu- 
lation at different redshifts, for different SMBH masses, 
and for a run including an approximation of optically 
thick cooling. The mass accretion rate is not steady, but 
fluctuates randomly and substantially, the accretion rate 
through the circumnuclear disk being almost as often 
ne gative as positive. This res ult is consistent with that 
of lHopkins fc Quataertl ()2009f ). who find a large amount 
of variation in the instantaneous accretion rates in their 
galaxy simulations on scales < 1 kpc. We find that there 
is no preferred timescale for accretion over the course 
of the ^million years spanned by each of the simula- 
tions. The flat slope of the Fourier transform of the ac- 
cretion rate characterizes the fluctuations as a function 
of timescale, revealing the stochastic nature of accretion 
in the simulations. This complex and chaotic behavior is 
apparent at each of the three redshifts explored in detail 
here: z = 3, 4, and 6, as well as in z = 4 simulations 
containing black holes with larger masses, and the z = 4 
run including optically thick cooling. 

The dynamic and turbulent nature of the simulated 
circumnuclear disk is not well modeled by Bondi accre- 
tion, even when using a modified prescription that takes 
into account the vorticity and turbulent properties of the 
gas. Bondi prescriptions do not take into account the ef- 
fect of the disk's self-gravity. Additionally, the modified 
Bondi accretion rates, estimated from the properties of 
the disk in the highest-resolution part of the simulation 
(where the density rises steeply with decreasing radius), 
are highly super-Eddington on the smallest scales (as- 
suming a radiative efficiency of 0.1). As the Eddington 
rate provides an upper limit for the accretion rate in the 
presence of radiative feedback from the SMBH, highly 
super-Eddington rates are not likely to endure in nature. 

The Eddington rate (or a few times the Eddington 
rate) is set as an upper limit to the accretion rate in 
many of the simulations that em ploy Bondi type pre- 
scriptions to mod el accr etion (e.g. JSpringel et"aLlf2 005a; 
Pi Matteo etaLl 120051: iLi et alj [20071: iDi Matteo et al.l 
20081 ) . On scales corresponding to the spatial resolution 
of those simulations (> 100 pc), the accretion rates pre- 
dicted by the modified Bondi prescription in the present 
simulations are closer to the Eddington limit than they 
are at smaller scales. However, black holes in other sim- 
ulations spend the majority of time accreting well below 
this limit, at substantially sub-Eddington rates. The dif- 
ference between the rates determined here and those of 
other simulations is even larger when comparing with 
prescriptions that multiply the Bondi accretion rate by 
a large factor to approximate the effects of the small 
scale physics, since no such factor has been included 
in our calculations (see discussion of the parameter a 
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iBooth fe Schavi[200l IJohansson et al.l l2009alTbh . How- 
ever, it is not straightforward to make a direct compar- 
ison between accretion rates in the present simulations 
and the rates in other simulations that contain prescrip- 
tions for AGN fee dback. As shown in galaxy merger 
simul ations (e.g., Springel ct al] l2005al DcBuh r et al.1 
2009), the presence of AGN feedback regulates black 
hole growth by forcing gas out of the circumnuclear re- 
gion. Future work incorporating AGN feedback into the 
present simulations will address this issue. 

The absence of AGN feedback in our simulations may 
result in an excessively dense circumnuclear gas disk, 
which corresponds to large, super-Eddington accretion 
rates within the Bondi prescription. Allowing the black 
hole particle to grow in situ in our simulations, or us- 
ing a more physically motivated recipe for star forma- 
tion in the zoom-in simulation would further change the 
distribution of mass in the circumnuclear gas disk. The 
circumnuclear disk is self-gravitating, almost certainly 
behaving o utside the Bondi regim e, even with the modi- 
fications of lKrumholz et all (|2006t ) . This result reinforces 
the need for high-resolution simulations incorporating 
detailed gas dynamics and radiative processes, in order 
to truly model the complicated dynamics on small scales, 
which ultimately govern accretion onto SMBHs. 

The chaotic behavior of the disk may be consis- 
tent with models suggesting that stochastic accre- 
tion of molecular clouds from the circumnuclear re- 
gion of the galaxy can p ower low-luminosity AGNs 
(Hopkin s fc Hernquistl[2006h . This manner of accretion 
may determine the spin of the black hole. We find that 
the angular momentum vector on scales < 100 pc can 
vary substantially from the direction of angular momen- 



tum on kiloparsec scales. In each simulation, the axis 
of the disk is closely aligned with the axis on kilopar- 
sec scales at the beginning of the zoom-in simulation, 
but slowly over time shifts direction as the circumnuclear 
disk dynamically evolves. 

The angular momentum can be misaligned with the 
rest of the circumnuclear disk, resembling the stochastic 
fueling scenarios mentioned in Section [5j As is observed 
in several of the simulations presented here, clumps of 
gas develop in the circumnuclear disk, which change the 
angular momentum of the gas that ultimately feeds the 
accretion disk. The change can occur as clumps fall into 
the center or as larger clumps develop within the disk 
and change the potential. However, as cautioned above, 
the simulation does not yet include feedback from an 
accreting black hole, which will heat the gas and drive out 
some of the material in the circumnuclear disk. Future 
work including more realistic physics in the simulations 
will allow us to constrain better the effects of accreting 
gas on the properties of the SMBH and the feedback it 
produces. 
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